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We present practical methods to measure entanglement for quantum simulators that can he realized with 
trapped ions, cold atoms, and superconducting qubits. Focussing on long- and short-range Ising-type Hamil¬ 
tonians, we introduce schemes that are applicable under realistic experimental conditions including mixedness 
due to, e.g., noise or temperature. In particular, we identify a single observable whose expectation value serves 
as a lower bound to entanglement and which may be obtained by a simple quantum circuit. As such circuits 
are not (yet) available for every platform, we investigate the performance of routinely measured observables as 
quantitative entanglement witnesses. Possible applications include experimental studies of entanglement scaling 
in critical systems and the reliable benchmarking of quantum simulators. 


I. INTRODUCTION 

Harnessing the potential of well-controlled experimental 
platforms, quantum simulators have recently emerged as ana¬ 
logue devices to study paradigmatic condensed-matter models 
[1]. To date, a considerable variety of devices have been pro¬ 
posed and partially realised to serve the central aim in this 
held, the preparation and control of quantum states with a 
number of constituents that is beyond the reach of classical 
simulations [2]. For the demonstration of genuinely quantum 
features of these simulators, it is thus of considerable interest 
to hnd methods which quantify entanglement and, if possi¬ 
ble, relate the hndings to classical simulatability. For pure 
states, the bi-partite block entanglement is a direct hgure of 
merit for the resources required when simulating many-body 
systems with numerical methods such as the density-matrix 
renormalization group [3-5]. One way to obtain the entangle¬ 
ment contained in a state in the laboratory would be to per¬ 
form full quantum state tomography [6] and to compute the 
entanglement of the reconstructed state. However, this is not 
only impractical due to the exponential resources required — 
the proverbial curse of dimensionality — but for many recon¬ 
struction schemes it may also lead to a systematic overestima¬ 
tion of the true entanglement content [7]. An experimentally 
feasible and rigorous alternative is to instead rely on lower 
bounds which may be obtained directly from measured ob¬ 
servables [25-30] and such lower bounds to the entanglement 
should (i) rely only on a few observables in order to avoid the 
curse of dimensionality, (ii) avoid assumptions on the state 
in the laboratory (such as, e.g., symmetries, temperature or 
an underlying Hamiltonian), and (iii) should be applicable to 
the experimentally relevant setting of mixed states. Indeed, as 
has already been demonstrated, (i)-(iii) may be met and entan¬ 
glement may be quantihed from signihcantly less observables 
than are required for the knowledge of the full state; E.g., col¬ 
lective observables are capable to detect [8-10] and quantify 
[11-14] entanglement. Note that extending (ii) also to observ¬ 
ables, is known as device-independent entanglement quantih- 
cation, for which there does not even need to be a quantum de¬ 
scription of the employed measurement device, see [15] and 
references therein. Here, however, we will assume that the 
relevant observables are actually those that are measured. 

We construct and analyze lower bounds to the bi-partite en¬ 
tanglement of states arising in the quantum simulation of a 


variety of spin models such as 

N N 

H=J2j.,,dlai+BY,ai, ( 1 ) 

i,j — l i—1 

which have recently been implemented in experiments with 
trapped ions [16-20], superconducting qubits [21], and ultra¬ 
cold atoms [22, 23]. We will consider ground states and 
their quasi-adiabatic dynamical preparation employing realis¬ 
tic noise models, including decoherence-induced mixedness. 

Our aim is to quantify bi-partite block entanglement of one 
part of the chain vs. the rest relying only on measurements 
of certain observables Ci. Denoting experimentally obtained 
expectation values of these observables by Ci [24], we are thus 
interested in 

£’mm[{C'z}, {Cj}] = mm |£;(p) I tr[(7ip] = Cij , (2) 

i.e., we consider the minimal amount of entanglement that is 
consistent with the obtained measurements c^. Here, E is the 
entanglement measure of choice and the minimization is taken 
over all density matrices g. As such, we follow the programme 
initiated in Refs. [25-30]. Note that no assumption on the 
state in the laboratory enters our considerations. While we 
will present tailored lower bounds to Smin that work particu¬ 
larly well—in some cases even providing E^nin exactly—for 
certain classes of states, we stress that all bounds presented in 
this work are valid for arbitrary states - pure or mixed. 

For systems governed by Hamiltonians as in Eq. (1), we 
identify a single key quantity in order to obtain lower bounds 
on i^min. That is, it turns out that a single observable C con¬ 
stitutes a common quantitative witness and, in fact, for large 
classes of states determines not only a lower bound but the en¬ 
tanglement of Q itself. We show how this witness may be mea¬ 
sured directly by employing a simple quantum circuit. If such 
a circuit is available, entanglement may thus be quantihed for 
systems consisting of an arbitrary number of spins. If it is 
not available, the above observation still allows us to trans¬ 
form the numerical minimization in Eq. (2) into the problem 
of computing the smallest eigenvalue of a sparse matrix and 
thus obtain results for more than 20 spins (and in principle 
many more using DMRG methods [5]). With recent imple¬ 
mentations of models as in Eq. (1) in mind, we thus introduce 
schemes for practical and rigorous experimental entanglement 
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estimation using only a few readily available observables and 
without relying on any assumptions on the state in laboratory. 

Throughout, we will use the logarithmic negativity [31] as 
our bipartite entanglement measure and consider the biparti¬ 
tion {1,..., ^}|{y + 1, • ■ • assuming N to be even. 
The logarithmic negativity is a full entanglement monotone 
for mixed states [32], an upper bound to the distillable entan¬ 
glement [33], and has an operational interpretation [34]. It 
reduces to the Renyi entanglement entropy with Renyi index 
1/2 on pure states, which, e.g., distinguishes topologically 
ordered phases (as do all the Renyi entanglement entropies 
[35]). In a setting involving mixed states, a topological con¬ 
tribution to the logarithmic negativity of the toric code model 
has been established in [36]. 


II. PRELIMINARIES 


We start by introducing the relevant quantities. The loga¬ 
rithmic negativity is defined as 

^^l.n.(p) = logllp^lll, (3) 

where is the partial transpose of g with respect to the cho¬ 
sen bipartition A\B (here, {1,..., -|- 1,... ,iV}) and 

||X||i = tr |X| = max{tr(C'X) | — 1 < C < l}. (4) 

is the trace norm. By its variational form we have that for any 
observable with — 1 < < 1 

Einig) > log(C)^. (5) 


Any observable C with this property thus serves as a quan¬ 
titative entanglement witness as it not only witnesses entan¬ 
glement but indeed provides a lower bound. As an impor¬ 
tant example for such a quantitative witness consider the un¬ 
normalized maximally entangled state 
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which fulhls —1 < ((7 |$) ($| < 1 for any unitary U = 

V 0 W. Hence, for any state g 

Ein.(g)>log max (C/|$)(T>| f/'l') (7) 

u^v^w 

The signihcance of the quantitative witness (7 |<i>)($| (7^ 
becomes clear when considering pure states; For a given 
pure state, consider its Schmidt decomposition = 

■*/'s|as)|fes) and let (7 = V 0W he the unitary that takes 

1^) to EJas)l&s)- Then = ||(|?/)(^/>|)f||y^ and 

thus Eq. (7) becomes an equality. 

While in general this requires the knowledge of the Schmidt 
vectors, we will see below that for large classes of states, 
equality may be achieved for one particularly simple uni¬ 
tary. This fact may be used to greatly simplify the opti¬ 
mization in Eq. (2). Eurthermore, for these classes of states. 


{U |‘i>)(4>| W)g may be obtained directly by applying a sim¬ 
ple quantum circuit as in Eig. 1 consisting of mutually com¬ 
muting N/2 two-qubit controlled-not and N/2 single-qubit 
gates and subsequently performing a projective measurement 
of |0)(0|®^ in the computational basis. 


III. RESULTS 


The Ising model in Eq. (1) has been realized on a variety of 
experimental platforms; Systems with tunable interactions are 
for example found in devices based on superconducting qubits 
[21]. Short-ranged couplings are encountered in experiments 
with ultra-cold atoms in optical lattices, see, e.g.. Ref. [22], 
in which nearest-neighbour interactions have been simulated. 
Eor ion-traps, the implementation of Eq. (1) has been pro¬ 
posed theoretically [37] and realised experimentally [16-18]. 
Here, the basic form of Jij is dictated by the properties of the 
trap and external laser fields; Eor the scheme demonstrated in 
[38], Ising couplings are generated by two non-copropagating 
laser beams with frequencies wq ± M where ujq denotes the 
energy splitting between the states defining a local spin, e.g. 
hyperfine clock states of ^^^Yb+ [18, 38, 39]. If the beatnote 
detuning fj, is sufficiently far from each (transversal) normal 
mode frequency uJm, in the usual rotating wave approximation 
and within the Lamb-Dicke regime this leads to an effective 
Ising Hamiltonian with couplings given by 
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where fli is the Rabi frequency of the ith ion, k the wave vec¬ 
tor difference of the laser beams, M the mass of the ions and 
hi m denotes the transformation between the vibrational site 
excitations and the normal modes [40]. The sum runs over all 
normal modes. The range of the interaction can be controlled 
by the detuning /i from infinite range if p is close to the cen¬ 
ter of mass mode frequency ujn, where all the spins couple 
equally to the motional degrees of freedom, to dipole-dipole 
interactions for p ^ ujn [39]. Alternatively the interaction 
range may be varied by changing the axial trap frequency ly^ 
[41]. In between these two regimes the couplings are well 
approximated by an algebraic decay. 


with 0 < p < 3. A transverse magnetic held may be intro¬ 
duced by an additional laser beam. Eurthermore, ferromag¬ 
netic couplings may be obtained by choosing different detun¬ 
ing fi or by initializing the system in the highest excited state 
and following an adiabatic protocol [18]. 


A. The Quantum Circuit 

Our main result is that for ground states of a variety of spin 
Hamiltonians the maximizing unitary in Eq. (7) may be given 
explicitly; 
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FIG. 1. Expectations of |^)(^| for the ground state of the ferromagnetic (left) and of for the ground state of the antiferro¬ 

magnetic (middle) long-range Ising Hamiltonian in Eq. (1) for realistic couplings [42] as in Eq. (8) with = 20 and B in units of 
Jo = 57 * \ Ji; i+i|/(A'^ — 1). The coupling range is determined by the detuning parameter /r. The expectations may be obtained via the 
circuit in Eq. (12), depicted on the right, and they coincide with the entanglement in the ground state, see Corollary 1. 


Corollary 1 Let H a Hamiltonian as in Eq. (1) and suppose 
it has a non-degenerate ground state |t/)). Let the couplings 
be such that Jij = If the "y x ^ matrix 

J with entries = Ji^Ar+i-y, i, j = 1..., negative 
semi-definite then 

Eln.{\f){'f\) = log 2 tr[|$)($| (10) 

If ff is positive semi-definite then 

Ein.{\i))ftp\) = log 2 tr[|$')($'| |V')(V'I]. (11) 

where |$') = ct® (g) • • • (g) (T ^/2 

This is a corollary of a theorem allowing for even larger 
classes of Hamiltonians which we prove in the appendix, 
where we also show that the conditions of the corollary are 
met by couplings as in Eq. (9). Furthermore, they are also 
met by the couplings as in Eq. (8) with, e.g., parameters as we 
choose them in the numerical examples. Hence, the bipartite 
entanglement (between the left and right half of the chain as 
quantified in terms of the logarithmic negativity) of any state 
that is a non-degenerate ground state of a Hamiltonian as in 
Eq. (1) with couplings fulfilling the hypotheses of corollary 1 
is equal to the expectation value of a simple (unnormalized) 
projector. What is more, this expectation value serves as a 
lower bound to the entanglement of any state (pure or mixed). 
One possibility to obtain this expectation value—so the over¬ 
lap of the state in the laboratory with |<i)), respectively |$')— 
is to apply a simple circuit and subsequently measuring the 
projector |0)(0|®'^: For the ferromagnetic case, (J < 0), we 
write |<i)) = R\0)^^, where 

N/2 

R = (^H^C,,N+l-^ (12) 

i=l 

and Ci j denotes the controlled-not gate acting on spin i (con¬ 
trol) and j (target) and Hi the Hadamard gate acting on spin 
i. The antiferromagnetic case, (J > 0), follows by addition¬ 
ally applying the transformation before the mea¬ 

surement. Note that in ion trap experiments, spin polarization 


measurements along a particular axis are routinely performed 
by spin-dependent resonance fluorescence. 

The logarithmic negativity of any state may thus be lower 
bounded by applying the circuit R, which is depicted in Fig. 1 . 
There, we also show numerical results for the thus obtained 
entanglement of the ground state of the Ising model in Eq. (1) 
for realistic [42] ferro- and antiferromagnetic couplings; cf. 
the phase diagram from the entanglement entropy in Ref. [43]. 

Let us emphasize again that R neither depends on the cou¬ 
plings Jij nor on the magnetic field B. Therefore, one does 
not require any knowledge about these parameters and the 
method is robust against an inexact implementation of the 
Hamiltonian as long as the hypotheses of Corollary 1 hold. 
Also, as |$)($| is a quantitative entanglement witness, the 
state in the laboratory g does not need to be exactly in the 
ground state, it does not even need to be pure, in order for 
to be a lower bound. 


B. Other observables 

Although experimentally feasible (see e.g. Ref. [44] for the 
realization of a CNOT gate in ion traps and Ref. [45] for super¬ 
conducting qubits), other observables may be more accessible 
than the implementation of the circuit R. To this end, we give 
lower bounds to Eq. (2) in terms of arbitrary observables Ci. 
Combining Eqs. (2) and (7), we find that i7min[{C'i}, {ci}] is 
lower bounded by the logarithm of the solution to the semidef- 
inite program (SDP) 

max > WiCi 

* , (13) 

subject to ''^^WiCi < |‘h)($| . 

i 

Considering this SDP instead of the original Eq. (2) leads to 
a significant simplification of the optimization problem and 
standard SDP solvers like, e.g., SeDuMi [46] may be used. 
Furthermore, the simplified SDP in Eq. (13) is directly acces- 
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sible to algorithms such as SDPNAL [47] or SDPAD [48] in¬ 
tended for solving large-scale SDPs with (real) matrices of di¬ 
mension more than 4000 and number of constraints of the or¬ 
der of 10®. Therefore these algorithms may outperform stan¬ 
dard interior point methods where they become too expensive 
computationally. Note that the observables Ci in Eq. (13) are 
entirely arbitrary and this scheme is thus sufficiently versatile 
to accommodate measurements of any experimental platform. 

Motivated by the fact that if the ground state is separated 
from the first excited state by an energy gap, the Hamiltonian 
itself provides an entanglement witness [49], we consider wit¬ 
nesses of the form 

N 

W = wol+(^al + wiH, (14) 

i=l 

where we included the (optional, see Fig. 2) operator o'* 
to account for the small gap in the symmetry-broken phase. 
This further simplifies the optimization in Eq. (13) as now we 
are considering only one observable—namely W —and the 
number of optimization variables is reduced to one. Note that 
for H as in Eq. (1), the witness W consists of at most quadrat- 
ically many observables (plus the single optional observable 
(tJ.) and hence its expectation value may in this sense 
be obtained efficiently: The experimental effort is reduced to 
obtaining the expectation value of the magnetization 
and all pairs it® (T^ for which Jij is non-zero. In ion-trap and 
superconducting-qubit experiments such observables are rou¬ 
tinely measured. For nearest-neighbour couplings as, e.g., the 
ultra-cold atoms experiment in Ref. [22], this amounts to only 
linearly many observables, the correlators may be 

obtained directly under a quantum-gas microscope [22], and 
the magnetization by a Fourier-transformation of the time-of- 
flight distribution. For the couplings one could either choose 
a theoretical prediction (for ion traps given in Eq. (8)) or, if 
possible, measure them experimentally (see the methods used 
in Ref. [19] for ion traps). Then the SDP may be avoided 
completely by choosing wq as the smallest eigenvalue of 

N 

( 15 ) 

i=l 

as then the constraint W < jT*) (<1>| is automatically fulfilled. 
As this operator is a sparse matrix, standard eigenvalue solvers 
allow for system sizes of more than 20 qubits. In fact, since 
|$)($| possesses a representation as a matrix product opera¬ 
tor of bond dimension four, DMRG algorithms may be used 
to obtain the smallest eigenvalue for much larger systems. In 
Fig. 2 we show numerical results for the above procedure. 
Again, we do not put any assumptions on the state in the 
laboratory—the expectation {W)g is a lower bound to the en¬ 
tanglement of any state g but, of course, we know that the 
bound will work particularly well for states that fall within 
the framework of Corollary 1, i.e., ground states of Hamilto¬ 
nians as in Eq. (1) with couplings as in Eq. (9) or as in Eq. (8) 
with parameters as for all the numerical examples considered 
here. 



B [J„] 


FIG. 2. Lower bounds to the entanglement of ground states of the 
Hamiltonian in Eq. (1) with A = 16, g — 117.6 kHz, Jo as in Fig. 1, 
and realistic couplings [42] as in Eq. (8). The hlack line shows the ex¬ 
act logarithmic negativity of the ground state. Entanglement bounds 
are obtained by optimizing the quantitative witness in Eq. (14) over 
wi. The couplings in the witness are as in Eq. (8) hut randomly per¬ 
turbed by 2% to mimic imprecise knowledge and shown are several 
random trials as density with the optional operator in blue 

and without in cyan. 


C. Quasi-adiabatic preparation and benchmarking 

In non-equilibrium situations, quantum simulators of one¬ 
dimensional spin systems may outperform classical comput¬ 
ers already for a moderate size of spins: As opposed to states 
in equilibrium, which typically have little entanglement (cf., 
area laws for ground and thermal states [3, 52]), the en¬ 
tanglement generated in non-equilibrium situations may be¬ 
come large [53]. Arguably the best numerical algorithms for 
the simulation of one-dimensional (non-)equilibrium quantum 
many-body systems are those based on matrix product states 
(MRS) and matrix product operators (MPO) [5, 54]. The re¬ 
sources required to treat such states numerically are directly 
related to their so-called bond dimension. For pure states, 
i.e. MRS, there is an intimate relation between the bond di¬ 
mension and the entanglement content as quantified in terms 
of Renyi entanglement entropies [4]. For mixed states, i.e. 
MPO, this connection is far less clear. Indeed, an MPO may 
have a small bond dimension while at the same time have 
a large block entropy—the product operator (1/2)®^ being 
the most striking example. In this sense, using pure-state en¬ 
tanglement measures (such as Renyi entanglement entropies) 
as benchmarks may lead to false conclusions because in ex¬ 
periments mixedness is unavoidable. We illustrate these rela¬ 
tions by considering the quasi-adiabatic preparation of ground 
states of Ising Hamiltonians as commonly performed in ion- 
trap experiments [18]: Initializing the system in a product 
state with all spins aligned parallel to the magnetic field, the 
field is reduced slowly (compared to the Ising interactions) 
until the desired B is reached. In a realistic setting, such a 
protocol is prone to noise processes such as non-adiabaticity, 
spontaneous emission (se) and dephasing (dph), which are 
considered the main noise sources [18]. We model this by the 
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FIG. 3. Quasi-adiabatic ramping of the magnetic field across the phase transition simulated using the model in Eq. (16) with Hamiltonian as 
in Eqs. (1) and (8) for A'’ = 8 spins and parameters as in [56], leading to a relatively long-ranged interaction ~ |i — and Jo = 2 tt ■ 3.3 

kHz. The magnetic field is ramped according to B = 1.1 • Jo, where k = 27r • 20 kHz and to = 0.6 ms. Shown are the block 

entropy ^(tri . . jv/ 2 [^(f)]) (an entanglement measure if the state was pure) in red, the exact bi-partite entanglement of the simulated state g{t) 
in black, lower bounds as obtained by the circuit R (blue), by the SDP in Eq. (13) with observables i,j = 1,..., N, a = x,y, z, 

and d^ • • • d^ as input (solid cyan), and by optimizing the quantitative witness in Eq. (14) over wi (dashed cyan). Shown in grey are upper 
and lower bounds [55] to the MPO approximation error in Eq. (17) (B = 1, 2, 3, 4 top to bottom). Note the monotonicity of the block entropy 
as opposed to the behaviour of the approximation error and the entanglement as quantified in terms of the logarithmic negativity. 


commonly used Lindblad quantum master equation 

^ + E 

i,<y 

■ 

with a = se, dph and Lf = , and 

{•, •} the anticommutator. Numerical results are summarized 
in Fig. 3 and the main conclusions are: The block entropy 
5'(tri..., Ar/ 2 [e(f)]) (a measure of entanglement if the state was 
pure) increases with time for all noise-strengths while the true 
entanglement reaches a maximum after which it decreases in 
time. From the block entropy one would thus falsely conclude 
that the state becomes harder and harder to simulate while the 
error when approximating g{t) by an MPO qd with bond di¬ 
mension D, 

eoit) = min\\g{t) - qdWf, (17) 

Qd 

reaches a maximum and then decreases in time [55] as does 
the entanglement. The exact mathematical connection be¬ 
tween approximability by MPOs, entanglement, and other 
quantities such as, e.g., mutual information, remains an open 
question however. 

IV. SUMMARY AND OUTLOOK 

In the setting of quantum simulations of the transverse- 
field Ising model, we have introduced methods to estimate 
bi-partite block entanglement without putting any assump¬ 
tions on the state in the laboratory. The principles pre¬ 
sented here are applicable to, e.g., ion-trap, cold-gases, and 


superconducting-qubit implementations and we have focused 
on the ion-trap platform for specific examples. A lower bound 
to the entanglement is given by the overlap with a certain state, 
which may, e.g., be obtained by a simple quantum circuit and, 
for large classes of states, actually gives the entanglement ex¬ 
actly instead of just bounding it. As obtaining this overlap 
may, depending on the platform, may represent a consider¬ 
able experimental challenge, we further investigated the per¬ 
formance of routinely performed measurements as means to 
estimate the entanglement. As we consider the benchmarking 
of quantum simulators as one possible application, we have 
compared the matrix-product-operator bond dimension, block 
entanglement, and block entropy for a quasi-adiabatic proto¬ 
col preparing ground states of the transverse-field Ising model. 
It is our hope that this inspires work towards revealing the 
exact mathematical connection between the MPO approxima¬ 
tion error and correlation measures and also towards schemes 
to estimate the former directly by measurements in a way sim¬ 
ilar to the schemes presented here for entanglement. 
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Appendix A 

In this section, we prove a more general version of Corollary 1. A central role is played by symmetry properties of the 
Hamiltonian. Notably, Hamiltonian (1) with, e.g., couplings in Eq. (8) or Eq. (9) is invariant if we either consider the spins 
from left to right or from right to left. More formally, we denote accordingly by I the transformation interchanging spins 
i -(A- N + l — i. Eirst however we consider SWAP invariant Hamiltonians, i.e. invariant under interchanging subsystems A ■(A- B. 
Below Theorem 1 we give the equivalent statement for models with I invariance. 

Before we are in the position to state the main theorem, we need the following simple fact; On a bipartite system Ha C) Hb = 
let (5 be a SWAP-invariant Hamiltonian, i.e. Sav^bGSa^b = G, then there are S R and operators A and Bi, such 

that 


(7=l(g)A + l(g)l+^ b^Bi (g) Bi. 

i 


(Al) 


To see this, note that any Hamiltonian can be written as 

G = 9i,jGi 0 Gj 

i,3 

= 1 (g) go.iGi + QifiGi (g) 1 + dijGi (g Gj, 

i i *5^^! 


(A2) 


with {1, (5i, (72) ■ • ■ } a hermitian operator basis and pij G R. If (7 is SWAP-invariant, we have po,i = 5i,o and the matrix g 
with entries j, j > 1, is real symmetric. Hence, orthogonal diagonalization OgO^ =: diag(6i) together with the definitions 

A:=J29oA, (A3) 

i 

Bj 'y ^ OijGi (A4) 

j 


yields the desired form (Al). 

If G is as in (Al) with A and Bi real (but not necessarily hermitian) and h<o for all i, we call G negative SWAP-invariant. 
Eurther, we let 


1 ^) ■= h"GN/2)\ii---iN/2) =:'^\i,i)- 


(A5) 


Theorem 1 Let G be a negative SWAP-invariant Hamiltonian with a non-degenerate ground state \^jt). Then ground state 
logarithmic negativity with respect to left half vs. right half is given by 


= logtr[|$)($| |^/>)(V^|]. 


(A6) 



8 


Proof. Since G is real and the ground state non-degenerate, we may write (up to a global phase) 

IV”) = I*’-?') (A7) 

id 

with 'ijjij real. Furthermore, since G is SWAP invariant we have Sa<^b |^) = IV') such that 

1^’ *) = X! 

i.i 

i.e., the coefficent matrix is real symmetric. Hence there is an orthogonal transformation O diagonalizing ip with diag(Ai) := 
OipO"'". This allows us to write 


fc V I /Vi / k 

such that |||V')(V'niy^ = Z^i l^il (^li/^) = Zfc ^fc- Now, by (Al) and as b i < 0, we find for the ground-state energy 

{ip\G\ip) = XiXj (oil (flil G\aj) \aj) (AlO) 

id 

= 2^ A- {ai\A\aP) -f A^A^ {ai\Bk \ajf (All) 

i k id 

>2^A2(a,|i|a,)+^6fc^|A,||A,|(a,|4|a,Z (A12) 

i k id 

= d\Gd) (A13) 


with \ip) := 'YPi |Ai| loi) \ai). As we assumed that the ground state is unique, we hence have Xi = e''^|Ai| such that 

\{m\ = Y.\^^\ = \mw^\d\ (A14) 

i 

which completes the proof. q 

The equivalence between SWAP and /-invariance may now be exploited to obtain a result for spin Hamiltonians which are 
invariant under /; If H is I invariant then IbHIb is SWAP invariant such that if IbHIb is negative SWAP invariant and the 
ground state \ip) of H unique, we have 

-Bi.n.(|'0)('0l) = logtr[/B|$)($|/B \ip){ip\] = logtr[|$)(T>| \ip){ip\] (A15) 

with 1$) as in the main text. In particular, the Ising Hamiltonian with transverse field in Eq. (1) is /-invariant if the couplings 
fulfil Jij = Jjv+i_i Ar+i_j. Hence, Corollary 1 follows from Theorem 1. Besides the transverse field Ising model mentioned 
here, other Hamiltonians allow for similar conclusions. For example the ground state of the XY-model without magnetic field 
and arbitrary anisotropy is also determined by the expectation value of a single projector. 

For a given coupling matrix it remains to show that it is non-positive (or non-negative) in the sense of Theorem 1 . In general, 
due to the symmetry the coupling matrix is a Hankel matrix. In the following we show that this is the case for algebraically 
decaying couplings (see Eq. (9)): 




1 

I*-j|^’ 


where p > 0. To simplify notation we define the A x A .matrix J as 


'^id ’— *^i 


i,N-\-l—j 


\N+l-i-j\p (^(*) + ((7))P > 


(A16) 


(A17) 


with ((i) = A±i _ i . Define d to be the matrix with entries dij = (C(0 + C0))~^- We show that the couplings fulfil the 
condition > 0 by a result of entrywise matrix calculus. To this end, let f[A] denote the matrix obtained from A by applying 
/ entrywise, i.e. {f[A])ij := f{Aij). According to Theorem 1.4 in [57], J' given by (A17) is non-positive for all powers 
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p > 0, if log[d] > 0 on D+ = {x G (D'^| = 0}. The proof of this statement follows from the fact that there is a constant 

T > 0 such that log[(i] + tE > 0, where the matrix E is defined by Eij = 1. To see this, we first rewrite log[(i] = f[E — d] 
with 0 < 1 — dij < 1, where f{x) := — log (1 — x). Hence, log[(i] + tE > 0 iff / [i? — d] < tE, i.e. iff 

f[E-d]< f[{l-e-^)E]. (A18) 

This follows from two observations: (i) / is a Schur-monotone (S-monotone) function [58], i.e. for two real symmetric matrices 
A, B with entries aij, bij G (—1,1), 


0<A<B^f[A]<f[B], (A19) 

and (ii) d is positive definite and hence we may choose r large enough such that d > e~'^E. To proof (ii), rewrite d using 

1 


di i — 


C(*)+C0') Jo 




(A20) 


Thus, for X G (D^, 


x^dx = 




df > 0, 


(A21) 


with equality iff = 0 for f S [0,1], i.e. only if a: = 0 since all C,{i) are distinct. Hence d is positive definite. Thus, 

since log[d] coincides with log[d] + tE on 79+, the desired positive semi-definiteness of log[(i] on 79+ follows. 






